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Confirmation of the VeLLO L1148-IRS: Star Formation at very low 
(Column) Density 
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ABSTRACT 

We report the detection of a compact (~ 5"; about 1800 AU projected size) CO outflow from 
L1148-IRS. This confirms that this Spitzer source is physically associated with the nearby 
(« 325 pc) LI 148 dense core. Radiative transfer modeling suggests an internal luminosity 
of 0.08 to 0.13 Lq. This validates L1148-IRS as a Very Low Luminosity Object (VeLLO; 
L ^ 0.1 Lq). The LI 148 dense core has unusually low densities and column densities for a 
star-forming core. It is difficult to understand how LI 148-IRS might have formed under these 
conditions. Independent of the exact final mass of this VeLLO (which is likely < 0.24 Mq), 
LI 148-IRS and similar VeLLOs might hold some clues about the isolated formation of brown 
dwarfs. 

Key words: stars: formation; stars: low-mass, brown dwarfs; ISM: globules; ISM: individual: 
LI 148; ISM: jets and outflows; infrared: stars 



1 INTRODUCTION 

The sensitivity of the Spitzer Space Telescope allows for searches 
and studies of young stellar objects that are much less lumi- 
nous than those studied before. While surveys of star-forming re- 
gions using data from the Infrared Astronomical Satellite (IRAS) 
were sensitive to sources with luminosities > 0.1 Lq (d/140 pc) 2 
(where d is the distance: iMvers et alJI 19871 ). Spitzer can detect ob- 
jects with luminosities as low as 4 x 10~ 3 (rf/140 pc) 2 , i.e. has 
a sentitivity better by a factor 25 jPunham et al.ll2008l) . Spitzer has 
imaged a large number of isolated dense cores thought to be star- 
less. About 60 of them were observed as part of the Spitzer Legacy 
Project "From Mo lecular Cores to Planet Forming Disks", or c2d 
jEvans et all2003l) . 

These observations led to a series of discov eries of Very Low 
Luminosity Objects (VeLLOs; for definitions see Kauffmann et al.l 



for L 1521F-IRS feourke et alj l2006h and L328-IRS jLee et all 
l2009h . These papers established the initial VeLLO sample and the 
first overviews of their pr operties. Other well-s tudied VeLLO can- 
didates are Cha-MMS l l lBelloche et a"I]|2006l) and L673-7-IRS 
jPunham et alj l2010al) . See Table [2] for an up-to-date list. The 
combined c2d survey data for core s and clouds reveals 15 candi- 
date VeLLOs jPunham et~ai]|2008h . Not all of them are well stud- 
ied, though. The first VeLLO ever d iscovered was IRAM04191. It 
was identified by Andre et al. ( 1999) on the basis of millimeter ra- 
dio observations (see Dunham et alj|2006l for Spitzer results). 



I2005L Idi Francesco et alj|200l iDunham et alj|2008h , infrared stars 
of luminosity ^ 0.1 Lq which are embedded in dense cores. 
VeLLO colors and envelope properties are consistent with those of 
class or class I protostars jLada|[l987l : lAndre et al.lll993l) . young 
stars in their main accretion phase that are deeply embedded in the 
dens e core from which th ey form. 

lYoung et alj d2004h prese nted the first Spitzer-identified 
VeLLO, L1014-IRS (al so see: ICrapsi et alj l2005bl : iHuard et all 
l2006l : lBourke et alii 20051) . Further detailed studies were presented 
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The evolutionary status and fate of VeLLOs is presently unclear 
and likely to be diverse. There is probably g eneral consent that 
VeLLOs have a present-day mass < 0.1 Mq jYoung et al .1120041 : 
also see Section 16. lb . However, some VeLLOs drive prominent 
outflows, and thus seem to accrete significantly. These objects 
have been interpreted as very young (~ 10 4 yr) protostars of 
stellar fi nal mass (e . g., IRAM04191, which has a prominent 
outflow; lAndre et all 1 19991 : iBelloche et ail |2002|) . In contrast, 
VeLLOs devoid of such outflows might not accrete at a high 
rate. Some were proposed to be proto brown dwarfs forming 
in isolation (i.e., not as low-mass components of clusters and 
binaries; e.g., L1014-IRS, which h as a ~ 3" [~ 750 AU] outflow 
only d e tected with interferometers : | Young et alj2004l : ICrapsi et alj 
l2005bl : lBourkeetal.1l2005l : IHuard et aljkood) . This concept has 
been bolstered by recent discoveries of binary brown dwarfs 
which are too loosely bound to have been e xpelle d by multiple 
stellar systems (e.g., iRadigan et ail 120091 and iLuhman et alj 
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Figure 1. The discovery of L1148-IRS. The LI 148 region contains two 
dark dust features evident at optical wavelength (top panel). These are 
also detected in maps of dust continuum emission (bottom right panel). At 
24 fim wavelength, the source L1148-IRS (framed by a box in all panels) 
manifests as an unusually bright point source (bottom left panel). It resides 
near a relatively faint (9 mjy per 11" beam at 1.2 mm wavelength after 
smoothing to 20" resolution) dust continuum emission peak (bottom right 
panel). Solid and dotted lines indicate the relative orientation of the differ- 
ent maps. See Figure [2] for further information. Coordinates are given for 
J2000. 



120091 : also see \protect \vrule widthOpt\protect\ 
href { http : //www . vlmbinaries . org/ } { http : // 
www . vlmbinaries . org / } ). 

The second Spitzer VeLLO to be disc overed — the source dis- 
cussed in the present pap er, L1148-IRS dKauffmann et alj|2005l ; 
first noticed bvlKunlf 19981) — always maintained a peculiar status. 
As shown in Figure Q] this object was discovered as an unusually 
bright Spitzer MIPS source (compared to the field stars), SSTc2d 
J204056. 66+672304.9, which is closely associated with a 1.2 mm 
wavele ngth dust emission pea k revealed by the c2d MAMBO 
survey iKauffmann et alj|2008l) . It differs from other VeLLOs in 
that no other such obj e ct is as inconspicuou s in dust emission 
jKauffmann et alj|2005t 120081 : IWu et alj|2007t) . This implies that 
L1148-IRS resides in a natal dense core of un usually low mass 
and a verage density (~ 0.15 Mq, few 10 4 cm~ 3 : [Kauffmann et al .1 
l2008h . 

Unfortunate ly, initial single-dish searches by 

IKauffmann et alj J2005h di d not reveal any molecular out- 



flow. Thus. lKauffmann etafl d2005h could not entirely rule out that 
L1148-IRS is jus t a galaxy (probabilit y ~ 10 -4 , using 24 /xm 
source counts by IPapovich et al .1 1 20041) or background protostar 
accidentally aligned with the LI 148 dense core. Formally, this 



means that L1148-IRS does not meet all of the iDunham et al.l 
d2008h criteria for confirmed VeLLOs. Confirmation would require 
an outflow detection towards L1148-IRS at the systemic velocity 
of the dense core. This would prove that L1148-IRS (1) is a YSO 
which is (if) physically associated with the LI 148 region. Here, we 
present an outflow meeting these criteria. 

Our study is organized as follows. Section[2]introduces the data on 
which our study is based. We discuss the discovery and character- 
ization of the outflow in Section [3] The physical properties of the 
embedded YSO are derived in Section [4] Information on various 
VeLLO natal cores are summarized in Section [5] This provides a 
catalogue against which LI 148 can be compared. In Section [6] we 
turn to a discussion of LI 148-TRS in the context of other VeLLOs. 
Section[7]summarizes our study. 

In order to improve the readability of this paper, some com- 
plementary material is removed to two appendices. Appendix [A] 
provides details on the photometric data used in this study. Ap- 
pendix iBlcomprehensivelv presents data f or which a cursory anal - 
ysis was executed in previous studies (i.e., IKauffmann et al.ll2005h . 
and describes new data indirectly supporting the outflow analysis 
executed here. 



2 DATA & OBSERVATIONS 
2.1 Distance & Previous Studies 



IStraizvs et all d 19921) derive a distance of (325 ± 25) pc fo r the 
LI 148 region. This is smaller than the 450 pc est imated by iKunl 
dl998l) . We adopt the result bv lStraizvs et al.l ( 1 19921) . s ince they fo- 
cus on a s: 1° x 1° field centered on LI 148, while IKunl i ll 9981) 
study an ~ 2° x 10° region in which LI 148 lies only at the edge. 
This is consistent with a preliminary 2MASS-based distance esti- 
m ate of (301 ± 55) pc (G . Maheswar & C.W. Lee, in prep.; based 
on lMaheswar et al]20icl) . 

The LI 148 dense core was included in the iLee et ak 1 cm 

12001U2004I) surveys for inward motions. No signs of such mo- 
tions were found, but the observed li nes were weak an d the po- 
sition of L1148-IRS was not probed. iKirk et alj d2005h observed 
LI 148 with SCUBA at 450 fim and 850 fim wavelength, but 
their small target field (~ 2' diameter) did not include the VeLLO 
candidate considered here. The L1 1 48-IRS s o urce w as then pre- 
sented by IKauffmann et alj J2005h . IWu et alj d2007l) did observe 
this source with SHARC-II at 350 fim wavelength, but failed 
to detect it. These data, co mbined with the chara cteristics of the 
Spitzer-detected source, led IDunham et al. I d2008h to the conclu- 
sion that L1148-IRS "has a high likelihood of being an embed- 
ded low-luminosity protostar" (their 'group 3'). The s ource is also 
includ ed in the 1.8 cm radio continuum surveys by IScaife et alj 
d2010h . L 1148-IRS is detected at a signal-to-noise ratio ~ 3. As 
shown bv lScaife et all and repeated in Table[2] the implied intrinsic 
luminosity at 1.8 cm wavelength is comparable to what is found for 
other VeLLOs (i.e. Fi.g cm d , the product of observed flux density 

and distance squared); 

Most recently, I Kirk et all d2009l) used data from the Spitzer 
Gould Belt Legacy survey to search the entire LI 148 region for 
young stellar objects (YSOs). Just a single candidate YSO, i.e. 
our L1148-IRS, is f ound. This dearth of YSOs is confirmed by 
IDunham et alj d2008T ), though they find one further candidate with 
rather low YSO probability (i.e., in their 'group 6')- Like neighbor- 
ing clouds, LI 148 does thus f orm stars on l y at a very small rate 
or efficiency dKirk et alj200^) . lNutter et alj d2009ri took Akari data 
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at wavelength of 65 to 160|J.m towards a ~ 0?5 x 0?5 field con- 
taining LI 148-IRS. They detect and confirm LI 148-IRS; they do, 
however, focus on the extended matter in this region and uncover 
evidence for external heating of the clouds. 

Detailed assessments of the pr operties of the natal den se core of 
VeLLOs were first presented bv lKauffmann et al.l J2003) . Based on 
maps in lines of CCS (from the the Effelsberg lOOm-telescope) 
and N2H + (from the IRAM 30m-telescope), the dense core has 
line widths at half intensity < 0.25 km s' 1 . Non-thermal motions 
are thus sub-sonic, since H2 at 10 K temperature has a line width 
of 0.44 km s" 1 . Towards LI 148-IRS, N 2 H + lines yield a systemic 
velocity of 2.60 kms' 1 . Significant velocity gradients and jumps, 
as well as multiple velocit y components in C 18 Q , are observed 
near the YSO. The letter bv lKauffmann et all J2005h unfortunately 
lacked the space to present the data in detail. These are therefore 

for the first time provided in Appendix iBl 

Based data from the IRAM 30m-telescope, [Kauffmann et al.l 

d2005h found no indications for significant chemical evolution. 
Specifically, CO is depleted by a factor less than 10, and N2D 4 " 
is not detected. See Appendix IB4I for more details. However, there 
is evidence for contraction motions within LI 148: line asymmetries 
indicative of such motions are observed towards all three positions 
highlighted in Fig. [2] These new data are presented in Appendix 

m 

The Kauffmann et al. I J2008I) dust emission maps (obtained us- 
ing MAMBO at the IRAM 30m-telescope) reveal an unusually 
small mass reservoir for a star-forming core. Division of the peak 
H2 column density of 6 x 10 21 cm -2 by the filament minor and 
major axis of the parental core (1' x 5' angular size; Figs. [T] and [2] 
Appendix lB2t suggest mean H2 densities of (0.4-2.0) x 10 4 cm' 3 . 
Comparison with model density distributions with central den- 
sity plateau and dln(o )/dln(r) = —2 for large radii (Fig. 5 of 
iKauffmann et alj|2008l) give a central density < 3 x 10 4 cm -3 . 
For a source-centered aperture of 4200 AU radius, a mass of only 
(0.14 ± 0.02) Mq is derivecQ. Division of this aperture mass by 
the aperture area gives an aperture-averaged column density of 
(Af(H 2 ))42oo au = 4.8 x 10 21 cm' 2 . 

These values are unusually low, even for starless cores 
( IKauffmann et alj2008l) . LI 148-IRS thus resides in a natal core of 
unusually low density and column density. In this respect, LI 148 is 
similar to the VeLLO core L1014, which was deemed to be "star- 
less" previous to Spitzer observations. This makes it so interesting. 

The mass concentration appears to be centrally peaked: as ex- 
plained in Appendix lB2l cuts through the emission suggest density 
profiles dln(e)/dln(r) ~ - 2. LI 148-IRS actually resides T! 6 
(i.e., 2500 AU) offset from t he dust emission peak. W hile this off- 
set is potentially interesting, IKauffmann et alj d2008l) use Monte- 
Carlo experiments to show that this dislocation is only marginally 
significant given the uncertainties of their data: the beam of 11" 
size, the low signal-to-noise ratio, possible anomalous refraction, 
and potential pointing errors of order 3" have to be taken into ac- 
count. In combination, these effects can lead to apparent offsets 



1 The aperture radius of 4200 AU was denned bv lMotte & Andrei fepOlh . 
In Taurus, this size conveniently corresponds to a radius of 30". Also, 
iMotte & Andrei argue that this radius roughly contains the mass that can 
be accreted in about 10 5 yr of free-fall collapse of an initially hydrostatic 
sphere of ~ 0.2 kms" 1 sound speed. However, except for these general 
considerations, the radius is somewhat arbitrary. We adopt this size to com- 
pare with previous observations. 



between the VeLLO and dust emission peak position, even if such 
offsets do not exist. 



2.2 Spitzer Space Telescope 

LI 148 was observed by the Spitzer Space Telescope with the In- 
frared Array Camera (IRAC; AORKEYs 5158144) and the Multi- 
band Imaging Photometer for Spitze r (MIPS; AORKEY 9440768) 
as part of the c2d Legacy program fevans et alj|2003l . PID 139). 
Further Spitzer observations (AORKEYs 14608128, 14615808, 
18160384) were acquired by the "cores2deeper" (deep IRAC and 
MIPS integrations; PID 20386) and "Going Long" (first MIPS 
160 /im images of the region; PID 30384) programs. I nitial ver- 
sions of these maps were used in lKauffmann et all l l2005h . The im- 
ages presented here are derived from new mosaics combining all 
of the aforementioned data. Independent imaging by the Spitzer 
Go uld Belt L e gacy S urvey (Allen et al., PID 30574) is discussed 
by iKirketall d2009h . Spectra of LI 148-IRS have been acquired 
too (PID 30563; AORKEY 18401280), but are not discussed here. 
Details of the o bservations, post-pr ocessing and source extraction 
are presented bv lEvans et al .1 l l2007h . 



2.3 IRAM 30m-telescope 

The IRAM 30m-telecope near Granada (Spain) was used to search 
for an extended outflow from LI 148-IRS in the (J = 2-1) tran- 
sition of 12 CO. Data were taken on 2004 December 19 (single 
pixel facility receiver; system temperature ~ 380 K in T m b -scale, 
mapped area ~ 2' x 2') and 2005 July 19 to 21 (HERA multi- 
beam array; ~ 900 K, ~ 4' x 4'). Th e single-pixel maps were 
already used in IKauffmann eta?] J2005h . The pointing corrections 
in azimuth and elevation were below 3", and the focus corrections 
below 0.45 mm, as checked by regular pointing and focus checks 
at least once every 3 hours. The intensity scale was converted to the 
main beam brightness temperature scale according to the efficien- 
cies reported on the IRAM websitfl Figure [2] gives an overview 
of the layout of these observations. The data were calibrated and 
analyzed using IRAM's CLASS software packaged 

In parallel to these outflow searches, we used a variety of other 
tracers to probe the physical and chemical conditions in the dense 
gas. This includes lines of N 2 H+ (1-0) and N 2 D+ (1-0), CS (2- 
1) and (3-2), C ls O (2-1), HCO+ and H 13 CO+ (1-0), and HCN 
(1-0). All observations were obtained in frequency-switched mode 
using the single-pixel receivers and the autocorrelator backends. 
The N 2 H+ and C 18 lines were mapped on 2004 December 19 
in the same ~ 2' x 2' field for which 12 CO data were obtained. 
Our N 2 H+ andC ls O maps were already used in lKauffmann et ail 
The remaining lines were observed on 2005 July 19 to 21 in 
relatively deep integrations towards the three positions near LI 148- 
IRS sketched in Figuref2](chosen to probe LI 148-IRS, its immedi- 
ate environment, and the neighboring dust core). Focus and point- 
ing uncertainties are as previously quoted for the outflow search. 
More details are given in Table IbTI CLASS was used to calibrate 
and exploit the data. We removed low-order baselines and partially 
used Fourier (FFT) filtering to further improve the baseline quality. 



2 |http : / /www . iram. es/IRAMES/| 

3 http : / /www . iram . f r/IRAMFR/ GILDAsTl 
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Figure 2. Overview of observing strategy and results. As shown by the 
grey shading in the bottom main panel, the region contains two parallel fil- 
aments detected in dust continuum emission. They have masses of 8.6 Mq 
and 3.1 Mq, respectively (a black solid line marks the adopted boundary 
between the filaments). The position of L1148-IRS and its outflow are in- 
dicated by dots and lines, respectively. Boxes centered on L1148-IRS in- 
dicate the regions searched for outflows. Circles mark positions probed for 
dense core contraction motions (Section [2.1l Appendix lB5t . The top zoom 
panel presents details of the MAMBO dust emission map in grayscale and 
contours (starting at, and spaced by, 2 mjy per 11" beam after smooth- 
ing to the 20" beam marked by a circle). Markers indicate the position of 
LI 148-IRS, its outflow, and of two of the positions searched for contraction 
motions. Note that LI 148-IRS is offset from the nearby dust emission peak 
(by a marginally significant amount). Coordinates are given for J2000. 

2.4 KPNO 4m-telescope 

As part of a larger program to survey the c2d cores in the near- 
infrared, deep J, H, and K s observations of LI 148 were ob- 
tained during September of 2004 using FLAMINGOS (Florida 
Multi-object Imaging Near-IR Grism Observational Spectrometer; 
lElstonlll998h on the 4m-telescope at Kitt Peak National Observa- 
tory (KPNO). We use the data acquisition and reduction strategies 
that have been described in lHuardetalJ l l200r3) . 

2.5 Plateau de Bure Interferometer 

The Plateau de Bure Interferometer (PdBI) observed LI 148-IRS 
on 6 occasions during the summer of 2006 (project P057; observed 
on May 13, June 5 and 15, July 22, August 2, and September 
1). Between 4 and 5 antennas of the array were used in the 'D' 
configuration (i.e., the most compact one). The 3 mm and 1 mm 
band receivers were tuned to the (J = 1-0) and (2-1) transitions 



of 12 CO, respectively ( 115.271202 GHz and 230.538000 GHz; 
IWinnewisser et all 1 19971) . All 8 units of the correlator took data; 
for each receiver, 2 units were centered on the line to acquire spec- 
tra at 0.1 kms" 1 and 0.2 kms -1 resolution, and 2 further units 
were symmetrically displaced from the line to collect continuum 
data with a bandwidth of 320 MHz. The primary beam sizes were 
20" and 40". The minimum projected baseline in D configuration 
is ~ 20 m, and so structures larger ~ 16", respectively ~ 33", are 
not at all probed by the maps. LI 148-IRS was in the phase center 
of the array. 

3C454.3 was typically used for bandpass calibration, while 
1928+738 and 2037+511 served as phase calibrators. As usual for 
the PdBI, MWC349 data yielded the initial flux calibration (as- 
sumed to have flux densities of 1.12 Jy and 1.7 Jy in the 3 mm 
and 1 mm bands, respectively). Based on this, the fluxes of the 
other calibrators were fixed and iterated to obtain an optimum cal- 
ibration across all epochs. The flux calibration can be checked by 
imaging the phase calibrator 2037+511 in those epochs where it 
is not used to derive this calibration (i.e., is not held 'fixed' in 
the calibration routines). We find that the observations reproduce 
the epoch-averaged mean flux density with maximum deviations 
of 7% and 20% in the 3 mm and 1 mm band, respectively. 

The data were calibrated and imaged (including cleaning) us- 
ing the GILDAS software suite. Here, we use all epochs to con- 
struct maps in the 3 mm band, but only the best 4 periods (i.e., 
excluding July 22 and August 2) are taken in the 1 mm band. This 
yields beam sizes of 4" 5 x 3" 4 (position angle east of north of 82°) 
and 2" 5 x l'.'7 (p.a. of 65°) at 3 mm and 1 mm wavelength, respec- 
tively. To improve the signal-to-noise ratio, the spectra are sampled 
into bins of 0.3 kms -1 width. The continuum data are clipped in 
the spectral range to exclude CO emission from our target source. 
This yields noise levels in the 3 mm and 1 mm band spectral line 
data of 23 mjy beam -1 and 46 mjy beam -1 , respectively, while 
we find 0.39 mjy beam" 1 and 1.0 mjy beam -1 for the contin- 
uum data. 



3 AN OUTFLOW FROM L1148-IRS 

Previous analysis failed to establish that LI 148-IRS (0 is a YSO 
and (if) is physically associated with the LI 148 region. Here, we 
show that the source drives a weak small-scale CO outflow (< 5", 
about 1800 AU in projection). This confirms the hypotheses (f) and 
("■)■ 

3.1 A Compact CO Outflow 

3.1.1 Single-dish Observations 

We initially used the IRAM 30m-telescop e to search for an outflow 
in the 12 CO (2-1) line. This search failed jKauffmann et alj2005h . 
It is briefly summarized below. 

The single-pixel facility receiver was used to map an area of 
2' x 2' size with 30" spacing in frequency-switched mode (layout 
indicated in Figure^. The top panels of Figure[3]show an example 
spectrum from this search: emission from the dense core (as well 
as "telluric" emissiorQ from Earth's atmosphere) are detected, but 

4 The LSR velocity of the telescope can, e.g., be found from 
the calculator provided on the homepage of the Far Ultravio- 
let Spectroscopic Explorer (switched to topocentric coordinates), 

|http : / /f use ■ pha ■ jhu ■ edu/support /tools /vlsr . html| 
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no obvious outflow signal was found (neither in the form of line 
wings, nor in terms of clear spatial trends in emission). However, 
Fourier components, as well as baselines of 2 nd order, had to be 
removed during data reduction. Thus, these observations were not 
sensitive to outflows manifesting in wide and faint emission lines. 

Therefore, we later raster-mapped a field of about 4' x 4' size 
with 24" spacing using the HERA multi-beam array in position- 
switched mode (indicated in Figure [2}. These data have very flat 
baselines. Unfortunately, the reference position was not free of 
emission. Therefore, the data are severely contaminated within 
about ±2 kms" 1 from systemic velocity. After smoothing to 
0.5 kms" 1 velocity resolution, these spectra have a noise level of 
0.8 Jy beam" 1 in the main beam brightness scale. Still, no outflow 
emission was found. Also the mean spectrum averaged over all map 
pixels does not reveal any CO emission beyond ±2 km s" 1 veloc- 
ity offset (at a noise level of 54 mjy beam" 1 ). 

Figure [3] compares the single-dish observations with the 
interferometer-detected outflow described below. For this compar- 
ison, the PdBI spectra in the 12 CO (2-1) and (1-0) transitions are 
integrated over source-centered apertures of 3" and 10" radius, 
respectively. The data show that the single-dish observations are 
sensitive enough to detect the outflow. The outflow emission is, 
however, massively confused wit h the dense core em ission. This is 
similar to the situation in L1014 dBourke et alj2005l) . 



3.1.2 Interfere/metric Obser\>ations 

The non-detection with single-dish telescopes did not rule out the 
existence of a compact outflow fr om LI 148-IRS, such as the one 
seen in L1014 dBourke et alj2005h . We therefore also used IRAM's 
Plateau de Bure Interferometer (PdBI) to observe the (1-0) and (2- 
1) transitions of 12 CO. These observations reveal an outflow from 
LI 148-IRS. As we detail below, it manifests in two symmetric ve- 
locity components (Figures [3] and O that form a bipolar outflow 
system (Figure |4j- 

In the PdBI-derived data cubes, significant 12 CO emission is 
detected in the 1.2 to 4.2 kms" 1 velocity range. ("Significant" 
means that, in a given position-position-velocity cell, the emission 
of both 12 CO transitions exceeds the noise level by a factor ^ 4. 
The velocity channels are spaced by 0.3 kms -1 .) The N2H + - 
derived systemic velocity is 2.60 kms" 1 (Section [2. 11 Appendix 
|B3I >. Relative to this we obtain maximum projected velocity offsets 
of —1.4 to + 1.6 kms" 1 for the 12 CO outflow emission. We ap- 
proximate this as ±1.5 kms" 1 in the following. Correspondingly, 
we use velocity ranges of 0.9 to 2.4 km s" 1 (blue-shifted lobe) and 
3.0 to 4.5 kms -1 (red-shifted lobe; we exclude the central chan- 
nel at 2.7 kms -1 for clarity) to calculate the velocity-integrated 
emission of the lobes. 

Figure [4] presents integrated intensities in order to explore 
the structure of the lobes. These are more apparent in the 12 CO 
(1-0) data. In this line, the blue-shifted lobe is shorter than the 
red-shifted one. The extent of the latter lobe is 5'.'5 ± 0'.'3, or 
(1800 ± 100) AU in projection; to measure this, we find the most 
elongated contour (at a signal-to-noise level of 9.5 ± 1.0), drop 
to 50% of this intensity, measure this contour's maximum offset 
from LI 148-IRS, and subtract the beam radius linearly. Manual 



and using the "source" command within the AS- 
TRO program of the GILDAS software package 
(http: / /www, iram. f r/IRAMFR/GILDAS/> . The LSR velocity 
of telluric lines is the negative of the telescope's LSR velocity. 
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Figure 3. CO spectra towards LI 148-IRS. The panels show 12 CO(J = 2- 
1) data (top) and (J = 1-0) spectra (bottom). For the (2-1) transition, the 
solid line presents the frequency-switched IRAM 30m-telescope spectrum 
(11" beam) taken towards L1148-IRS. The signal near 7 kms" 1 velocity 
is of telluric origin and thus not relevant here. The shaded histograms por- 
tray the data acquired with the Plateau de Bure Interferometer. These spec- 
tra are averaged over source-centered apertures of 3" and 10" radius for 
the (2-1) and (1-0) lines, respectively. The systemic velocity derived from 
N2H + observations (Section ^. 1| Appendix |B3> is marked by a dashed line. 



fitting of the bright and elongated outflow emission with a straight 
line (Figure [4]l gives a most likely position angle of 107° ± 22°. 
(To calculate the uncertainty, we derive the position angles when 
displacing the fitted line by half a beam at 5" 5 offset from LI 148- 
IRS.) This fit is guided by the fact that LI 148-IRS resides right 
between the respective intensity peaks of the two lobes; this sug- 
gests an outflow direction similar to the position angle of the line 
connecting the peaks. Unfortunately, though, the emission farther 
away from the star (> 5") shows no clear structure. This leaves the 
option of a wide-angle outflow propagating in a direction signifi- 
cantly different from the one obtained above. 

Figure [3] shows spectra of the outflow emission. Gaussian fits 
to the spectra, and comparison to the N2H + systemic velocity, 
yield projected velocity offsets of —1.0 and +1.2 kms" 1 for the 
12 CO (1-0) line, respectively -0.9 and +1.2 kms" 1 for the (2-1) 
line. We approximate this as ±1.0 kms" 1 in the following. Gaus- 
sian line fitting of the PdBI spectra gives velocity-integrated flux 
densities of 21.0 Jy kms" 1 and 8.9 Jykms" 1 (15% of the sinsle- 
dish flux density) for the (1-0) and (2-1) lines, respectively. 

The position-velocity diagrams presented in Figure [5] show 
that the lobes are kinematically well-separated along the entire out- 
flow axis. Some faint outflow features can possibly be traced to 
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Figure 4. The outflow lobes for both CO transitions. The greyscale 
map presents the emission integrated over the 0.9 to 4.5 kms -1 veloc- 
ity interval. Solid contours outline the emission of the blue-shifted lobe at 
0.9 to 2.4 kms -1 velocity, while dotted contours delineate the red-shifted 
lobe at 3.0 to 4.5 kms . These contours start at, and are spaced by, the 
noise level times a factor ±3; they are dashed for areas with negative emis- 
sion. The central cross marks the position of LI 148— IRS. Ellipses give the 
beam sizes at half sensitivity. The arrow marks the location and orientation 
of the position-velocity cut presented in FigurefS] 



velocity offsets as large as ±2.5 kms - . The CO (2-1) panel 
reveals that the knot at +5" offset is detected in several velocity 
channels, and is therefore probably real. The elongated lobe seen 
in the (1-0) maps appears to form a counterpart of this knot. 

The CO emission cannot come from a circumstellar disk. Kep- 
lerian motion at 1 km s _1 speed and 1800 AU radius would require 
a central mass of 2 Mq . This significantly exceeds the estimated 
mass of the central star (Section [6.1 1 ) and of the dense core mate- 
rial immediately surrounding the YSO (Section fc.ll l. The emission 
must thus come from an outflow. 

There is a remote chance that the compact CO emission is 
not from an outflow, but is just an artifact of the interferometer- 
induced spatial filtering of the dense core emission. We cannot rule 
this option out. However, as we will show in Section [331 the close 
association of CO emission and L1148-IRS alone is sufficient to 
reliably demonstrate the protostellar nature of LI 148-IRS. 




-10 -5 5 10 

Position [arcsec] 



Figure 5. Position-velocity diagrams along the trajectory indicated in Fig- 
ure|4] The spatial origin coincides with LI 148-IRS, and the offset increases 
in the direction of the arrow drawn in Figure \4\ Contours start at, and are 
spaced by, the noise level times a factor ±2. Dashed lines are drawn for neg- 
ative signal, and solid contours are used otherwise. The horizontal dashed 
line marks the systemic velocity derived from N2H+ observations (Section 
[2~T1 AppendixIMl 



3.2 Outflow Mass and Kinematics 

To derive the outflow mass, we assume the 12 CO emission to be 
optically thin, adopt an H2-to- 12 CO abundance ratio of 10 4 , a 
mean mass p e r H2 m olecule of 2.4 H-masses, and use Eq. (A4) of 
iBourke et al. I i 19971) to calculate the column densities (for 12 CO, 
we use a permanent dipole moment of 0.11 D and a rotational con- 
stant of 57.636 GHz). We adopt an excitation temperature of 20 K. 
This is a bit more than canonical dense core temperatures ~ 10 K, 
but less than the few 10 2 K observed in outflo ws of much higher ve- 
locity (see, e.g. JSantiago-Garcfa et alj|2009l for a summary). This 
temperature uncertainty introduces a moderate mass uncertainty: 
temperatures in the range 10 to 50 K would give masses different 
by less than a factor 2. We obtain outflow mass estimates, m out , of 

1.3 x 10 -3 M Q and 3.9 x 10~ 4 Mq for the 12 CO (1-0) and (2- 
1) observations, respec tively. Correctio n for typical 12 CO outflow 
optical depths of 2 to 5 dLevreaultl 19881) would increase the masses 
by similar factors. For the (2-1) line, comparison with single-dish 
data implies a total flux loss of a factor 7 (Figure[3j. This could lead 
to underestimation of the outflow mass by similar factors. 

We adopt a characteristic outflow velocity of v ut = 
I ± 1 kms _1 |, and a radial size of r — 1800 AU. These are, of 
course, projected properties, which in principle must be corrected 
for the inclination angle. This angle is, unfortunately, not known, 
which introduce s systematic uncertaintie s int o our analysis . Usin g 
the formalism of lCabrit & Bertoutl dl990h and lBourke et all i 1997b . 
we derive the properties listed in Table Q] The quoted uncertain- 
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Table 1. Outflow properties. The quoted uncertainties refer to the range in 
outflow masses derived from the different 12 CO transitions. 



Property 


Definition 


Value 


size (radius) 


r 
AU 


1800 


velocity 


"out 

kms -1 


1 1 -1 1 

1 ± 1| 


mass 


m out 
10" 3 M 


n a 1 q 

U.4 tO l.o 


momentum 


"lout ■ V ou t 

10~ 3 Mq kms- 1 


0.4 to 1.3 


energy 


«W ■ -"out/ 2 
10- 4 Mq km 2 s~ 2 


0.2 to 0.6 


mechanical luminosity 


m ou t ■ vl u J[2r] 
10-5 L 


0.4 to 1.3 


force 


"lout • ^out / 1 " 


0.5 to 1.5 


10- 7 Af o kms- 1 yr- 1 


dynamical age 


r/v ut 
10 3 yr 


8.5 


rate 


m out ■ v out /r 
lO- 7 M yr- 1 


0.5 to 1.5 



ties refer to the range in outflow masses derived from the different 
12 CO transitions. Given their linear dependence on m ut, the out- 
flow properties also suffer from the systematic uncertainties affect- 
ing mass estimates. 

The properties of LI 148-IRS and its outflow are unusual, but not 
extreme for an embedded YSO. LI 148-IRS is l ess luminous than 
any source in the comprehensive lWu et alj j2004) outflow compila- 
tion. It also has an outflow mass, mechanical luminosity, and force 
smaller than observed for all (but 4, allowing for uncert ainties) of 
the 292 outflow systems with known luminosity in the IWu et al.l 
study. These trends do not merely reflect the low luminosity of 
LI 148-IRS; for LI 148-IRS, the above properties are still 1 to 2 
magnitudes below w hat is derived whe n extrapolating typical out- 
flow properties from IWu et al. I J2004) to the internal luminosity 
of LI 148-IRS. These deviations are, however, marginally within 
the scatter observed for other YSOs, and the LI 148-IRS outflow 
is thus not extreme. T he outflow survey of embedded YSOs by 
iBontemps et al .1 dl996t) can be used for another comparison. Based 
on this study, the force of the LI 148-IRS outflow is consistent with 
what is expected for class I sources of comparable internal luminos- 
ity. The YSO envelope mass of LI 148-IRS does, however, exceed 
the typical one expected for the observed outflow force by a factor 

~ 10. 

lAndreetalJ dl999j IRAM041 91), iBourke et alj feOOl 
L1014-IRS) and lPhan-Bao et al.l d2008l . ISO-Oph 102) report out- 
flow properties for other VeLLOs and brown dwarfs. The outflow of 
LI 148-IRS is similar to those of L1014-IRS and ISO-Oph 102 in 
terms of velocity (all have ~ 1 km s -1 ), but its mass 
vs. ~ 10~ 5 Mq, without mass corrections) and size (1800 AU vs. 
^ 1000 AU) appear to be larger. This is small, though, compared 
to the ~ 10 kms -1 velocity range and ~ 15000 AU extent of 
IRAM04191. 



3.3 Absence of Compact Dust Emission 

Compact dust emission, such as emitted by circumstellar disks, 
is not detected towards LI 148-IRS. At 1 mm and 3 mm 
wavelength, we observe intensities of 0.23 mjybeam -1 and 
0.56 mjybeam -1 , respectively. This is insignificant compared to 
the noise levels reported in Section [231 If we add the noise level 
times a factor 3, we obtain upper flux density limits of < 3.2 mjy 
and < 1.7 mjy, where we assume that the dust emission should be 
more compact than the beam. 

These upper limits are in line with, e.g., interferom etric con- 
tinuum observations of L1014— IRS dBourke et alJl2005T) . To derive 
mass l imits from these flux den sities, we adopt the framework laid 
out bv lKauffmann et all d2008h . Assuming dust with temperatures 
^ 10 K, this yields mass limits < 0.02 Mq and < 0.13 Mq, 
respectively. 

3.4 Near-Infrared Nebulosity 

The near-infrared (NIR) FLAMINGOS imaging observations pre- 
sented in Figure |6]reveal elongated nebulosity towards LI 148-IRS. 
This nebulosity is likely related to outflows from the star. 

LI 148-IRS is detected in all bands, but we only use the two 
bands with the best data (i.e., H and K s ) for our analysis. Fits of 
elongated Gaussians to the //-band image (which has the higher 
SNR) give a position angle (east of north) of 133° for the major 
axis. Along this axis, //-band emission from LI 148-IRS above 
10% of the peak intensity has a width of 4'.'l (3'.'7 in the K a band), 
while nearby stars in the image have a width of only 2" 5 (2'.' 7 for 
the K s band) for this position angle. Repetition of this procedure 
for the minor axis, and removal of the instrument-induced source 
width at 10% peak intensity, suggests an intrinsic extension of order 
3'.'2 x 2"4 or 1040 x 780 AU 2 . Both bands reveal an arc-like 
feature about 5" to the west of LI 148-IRS. This could be related 
to the aforementioned western blob in the red-shifted 12 CO (2-1) 
outflow lobe. 

The CO-detected lobes and the NIR nebulosity are well 
aligned (133° vs. 107° ± 22°). The nebulosity thus seems to trace 
structures such as an outflow-blown cavity. 

3.5 L1148-IRS as a Confirmed VeLLO 

As pointed out bv lDunham et alj d2008h . confirmation of a VeLLO 
candidate requires (i) that the internal luminosity is below 0.1 Lq, 
and {it) proof of emb eddedness in a dense core. In the character- 
ization framework of iDunham et alj {2008), the o uflow discovery 
furnishes the latter. Following the terminology of IDunham et all 
LI 148-IRS is thus a "confirmed embedded low-luminosity proto- 
star" (their 'group 1'). 

However, given the relatively weak emission, low velocity, 
and small lobe separation in LI 148-IRS, one may doubt that the 
compact CO emission is really coming from an outflow. As we 
show here, though, the proximity of compact CO emission and IR 
source independently demonstrates that LI 148-IRS is physically 
associated to the dense core. Then, LI 148-IRS must still be deeply 
embedded, since there is no other way to produce the Spitzer MIPS 
emission at the observed intensities (Section[4]l. 

We establish the association by computing the chance align- 
ment probability between background sources (i.e., galaxies) and 
the 12 CO emission. A background source density ~ 2000 sr 1 is 
to be expected for 24 /im sources brighter 70 mjy dPapovich et alj 
120041) . LI 148-IRS and the 12 CO (2-1) emission peaks are aligned 
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Figure 6. Images of L1148-IRS in the H (bottom panel) and K s bands 
(top panel) as seen by FLAMINGOS, shown as their negative. Contours 
delineate the 12 CO (1-0) outflow lobes for signals exceeding the noise by 
a factor 6 (see Figure |4j. The arrows indicate positions and directions of 
the cuts through the intensity distribution that are discussed in the text. The 
images reveal that LI 148-IRS is elongated, with the major axis at position 
angle RJ 133° . In addition there is an arc-like extended feature about 5" to 
the west of LI 148-IRS. To facilitate comparison between both images, the 
arc is indicated by a solid tracer line. 



to within 1" distance (Figure 0- For an individual cloud, the prob- 
ability for chance alignment is thus it (l") 2 • 2000 sr _1 = 1.5 x 
10~ 8 . For the entire c2d MAMBO survey, with ~ 110 dust emis- 
sion peaks, the probability of one chance alignment is 1.6 x 10~ 6 
(assuming that every dust emission peak contains interferometer- 
detected CO structures like those presented here). There is no re- 
alistic option that dense core and LI 148-IRS are not physically 
associated. 



4 THE SPECTRAL ENERGY DISTRIBUTION OF 
L1148-IRS 

The outflow observations establish the protostellar nature of 
LI 148-IRS. However, we still need to characterize the embedded 
infrared source. Here, we gather photometry of this object, sketch 
a model for the analysis of the observations, and discuss the prop- 
erties derived using this model. 
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Figure 7. The spectral energy distribution of LI 148-IRS. Bars represent 
flux density measurements with their associated uncertainty and filter width, 
while arrows indicate upper limits. Solid and dashed lines describe a dense 
core model containing a point source with 0.1 Lq internal luminosity and 
3500 K effective temperature; the solid line holds for an aperture of 2" 
radius, which is used to model the observations at ^ 10 urn wavelength, 
while the dashed line refers to an aperture of 4200 AU radius, as applied to 
the observations at larger wavelengths. This is detailed in Section l4~2l 



4.1 Observations 

Figure[7]presents the spectral energy distribution (SED) of LI 148- 
IRS. T hese data are assemble d from the 2MASS point source cat- 
alogue dSkrutskie et all2006l) . our FLAMINGOS observations, the 
IRAS faint source catalogue jMoshir&alJll99Ch . the final release 
of the c2d source catalogue l lEvans et alj|2007l) . our cores2deeper 
Spitzer project (PID 20386), our deep MIP S 160 Aim surve y of 
cores (PID 30384), SHARC-II i maging by IWu et all llOOH) , the 
SCUBA legacy source catalogue dFrancesco et al. 20081) . and our 
c2d MAMBO survey dKauffmann et al.ll2008l) . Appendix lAl dis- 
cusses how we consolidate these data from different origins. Table 
lAll lists the flux densities in detail. 

Integration of the SED yields bolometric luminosities and 
temperatures of 0.12 Lq and 110 K, respectively. We derive these 
parameters using power-laws to interpolate between the detected 
flux densities. We further ass ume flux densities scalin g with wave- 
length as A 4 for A > 350 /an jKauffmann et alj2008l) . Propagation 
of the flux-weighted observational uncertainties gives uncertainties 
< 20% in luminosity, and ±15 K in temperature . This is in agree- 
ment with t he LI 148-IRS characterization by Kauffman n et al .1 
d2005l 120081). bu t a bit warmer and fainter than what is found by 
iKirk et all J2009I : probably because they use MAMBO fluxes for 
larger apertures). 



4.2 Modeling 

We model the SED of LI 148-IRS as a spherical dense core that 
contains a point source and is externally illuminated by the inter- 
stellar radiation field (ISRF). The point source is assumed to emit 
like a black body of luminosity Lint and temperature T e g; we refer 
to these as the internal lumi nosity and the s tellar effective tempera- 
ture. The ISRF is taken from |Peraultl < ll987t) . We reduce its intensity 
at the core surface by a factor 4. This reflects shielding by extended 
material of a total (surface to surface) extinction of 3 mag, and is 
necessary to simultaneously match t he MAMBO and MIPS 3 obser- 
vations. We adopt dust opacities as IWeingartner & Draind l l200ll) 
predict for their "case B" and a total-to-selective extinction ratio 
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of Rv =5.5 (since this model provides a good fit to observed dus t 
opacities at wavelengths < 10 /im; e.g., Ilndebetouw et alj|2005h . 
Model flux densities are calcul ated using a code k indly provided by 
E. Kriigel (e.g., Appendix A of lPillai et af]|2006h . 

The dense core's density profile, g(r), is modeled as a con- 
tinuous step-wise power law. We truncate the sphere at an outer 
radius of 32500 AU « 0.16 pc, corresponding to « 100". The 
density is assumed to be negligible within the inner truncation ra- 
dius, rin. For most radii, the slope of the density profile is set to 
d ln(g)/d ln(r) = —2. However, good matches of the SED re- 
quire a flatter density law in some inner part of the sphere. We 
do this by introducing a constant density between the inner and 
outer transition radii, rj" ans and r°^ ns . Our fits use r°^ ns = 
2 x 10 14 m « 1340 AU (i.e., 8'.'25 diameter), since this is large 
enough to provide good SED fits (as demonstrated by the results 
below) but is compact enough to be consistent with the steep dust 
emission intensity profile that rules out extended regions of con- 
stant density (Section 12.11 ). With these assumptions, the density 
profile is fully defined once the density at the outer truncation ra- 
dius, gout, is specified. We fix the latter by requiring the model to 
exactly match the MAMBO flux density measurement. This even- 
tually leaves two envelope parameters (i.e., n n and rj" ans ) and two 
source parameters (i.e., Li n t and T e ff) unconstrained. 

At wavelengths < 10 /im, we attempt to minimize the mean 
squared relative difference between the observed and predicted flux 
densities (i.e., the difference between model and observation, di- 
vided by the observation; this property, x? e i> takes the role of the 
classical \ 2 )- We do not normalize the differences by the obser- 
vational flux density uncertainties (as usually done to derive x 2 ) 
because the predicted flux densities do sens itively depend on the 
adopted geometry (e.g., IWhitnev et ai1l2003l) . In this situation, the 
deviations between good models and observed flux densities are 
not dominated by observational uncertainties. We further require 
that a valid model does strictly not exceed the J-band limits. The 
synthetic observations are calculated for an aperture of 2" radius, 
and they are corrected for an intensity background determined be- 
tween 5" and 10" radius. This is used to capture the effect of the 
background subtraction in the 2MASS, FLAMINGOS, and IRAC 
observations. 

At wavelengths > 10 /im, we require the model to strictly 
match the observations within their observational uncertainties (ex- 
cept for all IRAS and the MIPS3 bands, which are too uncertain). 
This appears appropriate since in this wavelength domain, devia- 
tions between good model and observation should be dominated 
by observational uncertainties. We use an aperture of 4200 AU ra- 
dius to calculate synthetic observations. This is exactly the aperture 
used in MAMBO observations (Section |2.1t . The exact aperture 
size does not matter for the MIPS1 and MIPS2 bands, since the 
emission originates (in the model) at radii smaller 4200 AU. 

The PdBI-derived upper limits for dust continuum emission 
are not useful for our analysis of the SED of the LI 148 core on large 
spatial scales. Specifically, the interferometer probes radii ^ 4"5, 
i.e. ^ 730 AU, while the analysis above considers scales as large as 
4200 AU. The PdBI maps are only useful to examine the (absence 
of) disks around LI 148-IR, as done in Sec. 13.31 
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Figure 8. Quality of SED fits in dependence on the model effective temper- 
ature and internal luminosity. Only cells with shading or crosses have been 
examined. Grey shading gives X 2 c p the mean squared relative deviation be- 
tween SED model and observations at wavelengths < 10 (im. Crosses are 
drawn if no model can fit all of the SED constraints. As explained in the 
main text, only solutions with x 2 e i < 0-3 are deemed acceptable. These 
solutions are enclosed by a black line. The modeling thus constrains the 
internal luminosity to 0.08 < L lnt / 'Lq J$ 0.13 and the effective tempera- 
ture to 2500 < T cff /K < 5000. 

to a mean deviation of V0.08 = 28% between the observations 
and the model for wavelengths < 10 (jm. Here, we require that 
good fits not to exceed this mean deviation by more than a fac- 
tor 2, i.e., we require Xrel < 0-30 ~ (2 x 0.28) 2 . For certain 
luminosity-temperature pairs, no solution exists, since no model 
simultaneously fits the strict J-band and MIPS constraints. 
The modeling constrains the internal luminosity to 

Lint = 0.08 to 0.13 Lq. 

LI 148-IRS does thus marginally fulfill the VeLLO luminosity cri- 
terion of Li n t ^0.1 Lq. The luminosity is well constrained, since 
it only requires a good SED match and does not sensitively de- 
pend on model parameters. This is different for the effective photo- 
spheric temperature, which we find to be T e s = 2500 to 5000 K. 
The temperature basically follows from the SED slopes at wave- 
lengths < 10 (im. Unfortunately, at these wavelength the spherical 
geometry of our model is probably too simplistic to capture the 
relevant inner envelope extinction and scattering processes. This 
induces further systematic uncertainties. In any event, though, a 
central source of several 1000 K must be present to explain the 
detections in the 2MASS/FLAMINGOS bands. 

The SED of LI 148-IRS is thus consistent with the idea of 
LI 148-IRS being embedded in the dense core. In particular, this 
model provides a quantitative explanation for the failure to detect 
the VeLLO with SHARC-II. Initially, this lack of bright 350 jixm 
emission was interpreted as evidence against LI 148-IRS being 
an e mbedded VeLLO, since all other V eLLOs show this emis- 
sion dWu et al.ll2007l ; lDunham et al.l2008h . Our SED model shows, 
though, that no emission is to be expected, even if an embedded 
source is present. Slightly more sensitive 350 /im observations 
should be able to detect emission, though. 



4.3 Results 

Figure [8] presents the smallest mean squared relative deviation 
found (by varying n n and rj" ans ) for a given pair of luminosity and 
temperature. We find a minimum value of 0.08. This corresponds 



5 VELLO NATAL CORES AND OUTFLOWS 

Comprehensive understanding of LI 148-IRS and other VeLLOs 
requires a good knowledge of the environment in which these ob- 
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jects form. Therefore, we launched an investigation of the proper- 
ties of the LI 148 dense core. These results, as well as those for 
other VeLLOs, are collected here to support our further discussion. 
We also include additional data on outflows. Sections 15.31 and [5741 
summarize this information. 



5.1 Data for L1148-IRS 

Section |2T| summari zes the data situat i on for LI 1 48. The informa- 
tion is collected from l Kauffmann et al. I ( 120051) and l Kauffmann et al] 
d2008l) . As described in this summary, some of the data are also pre- 
sented in Appendix|B] Please refer to those sections for details. 



5.2 Data for other VeLLOs 

Here, we review the properties of known and proposed VeLLOs. 
First, we assess the quality of VeLLO candidates. Specifically, we 
consider objects that have previously been discussed to be similar 
to L1014-IRS, have indications for luminosities sj 0.1 Lq, or were 
discussed in the context of Brown Dwarf formation. 

Visual inspect ion of Spitzer images of Oph B-ll 
dGreaves et ai1l2003l) does not show any far-infrared counterpart. 
The core may actually be starless and is not co nsidered further 
here. We also tentativel y remove Lupus 3 MM S dTachihara et al.l 
12007b and L1415-IRS dStecklum et al] |2007|) from the VeLLO 
category. In these cases, estimation of the internal luminosity on 
basis of the flux density at 70 fim wavelength dDunham et al] 
12008b suggests luminosities significantl y exceeding the V eLLO 
limit of 0.1 Lq. We rate IC1396A:7i, dReach et al] 12009b as a 
source of unclear type, since the above luminosity estimate is 
not possible for this source. Preliminary photometry of the VLA 
source J02 5616+192703 in the L1457 dense core implies a low 
luminosity dHeithausen & B6ttnej201oh . but data are too limited to 
judge t he probability t hat th is really is a YSO and not a background 
object. iBarrado et al] d2009t) pr esent SSTB213 J041757, a p ossible 
VeLLO candidate in Taurus. iLuhman & Mamaie k (2010) show 
that this is unlikely to be a Taurus YSO, though, and SSTB213 
J041 757 is thus exclud ed here. 

IChen et al] d2010h and lEnoch et all d2010l) hav e recently pre- 
sented good candidates for "first hydrostatic cores" dLarso 5 1 1969b 
with internal luminosities < 0.01 Lq. These could be dense cores 
which have just become optically thick to their own radiation, be- 
fore eventually collapsing into a more compact YSO. Such ob- 
jects have a low age and luminosity, and they are embedded in a 
dense core . T his qualifies fi rst hydrostatic cores as VeLLOs. The 
IChen et all and lEnoch et al] sources will be included in our compar- 
ison once their nature has been established by decisive evidence. 

Table [2] lists 7 well-studied VeLLOs in order of decreasing 
aperture mass, AL^qqa u- It evaluates some criteria for dense core 
evolution forwarded bv lCrapsi etTaTJ d2005ah . By comparing prop- 
erties like column densities of dense gas tracers to characteristic 
thresholds, these criteria essentially test whether a dense core con- 
tains a significant mass of high density gas. Where appropriate (i.e., 
not when evaluating the CO depletion factor: Appendix IB4t . we 
convert published value s to match the dust properties adopted by 
lKauffmannetal.ld2008l : their Appendix A). We also list internal lu- 
minosities and results of outflow searches, as reported in the respec- 
tive discovery paper. Dista nce-normalized flu x densities at 1.8 cm 
wavelength, as reported bv lScaifeetaIld2010h , characterize the ra- 
dio emission properties of the VeLLOs. 



5.3 Properties of VeLLO Natal Cores 

The VeLLOs in Table [2] can roughly be divided into two 
groups. The first group is composed of LI 52 IF, Cha-MMSl, and 
IRAM04191. These cores fail in none of the lCrapsi et~ai1 d2005i]) 
criteria for "evolved" dense cores (if we stick only to criteria that 
can be tested with existing data). For instance, these cores are 
very dense (in terms of volume and column density), they have 
strongly peaked density profile (as indicated by small values of 
r 70 %, i .e., the radius of the dust emission contour at 70% peak in- 
tensitv; lcrapsi et al.l2005al) . and the inner regions of the dense core 
are in an advanced stage of chemical evolution (i.e ., CO de pletion, 
and production of deuterated N?H + ). ICrapsi et all d20053) argued 
that their criteria select dense cores relatively close to the onset of 
protostellar collapse. It is thus not surprising to find that collapse 
has indeed started in some of these objects. 

The second group is c omposed of " unevolved" dense cores, 
that fail in at least one of the lCrapsi et al] criteria; these cores have 
relatively low (column) densiti es, and there are ofte n no signs of 
advanced chemistry. Following ICrapsi et al] d2005j) , these dense 
cores are not expected to collapse and form stars. Section [63] ex- 
plores how stars might still have formed in these cores. The most 
prominent member of this group is L1014; the discovery of a 
Spitzer source in this core of lo w density, which w as believed to 
be "starless", was a big surprise dYoung et al]|2004l) . Also LI 148 
belongs to this group. 

The terminology used above is not meant to indicate that "un- 
evolved" cores mature and eventually b ecome "evolv ed" ones! It 
is only used to indicate the relation to the lCrapsi et al] criteria. 

In summary, some VeLLOs fulfil the lCrapsi et ail d2005al) criteria 
for evolved cores, while others do not. Since star formation is on- 
going in all cores listed in Table[2] the initial conditions for stars to 
form appear to be very diverse. 

L1148-IRS different from the other VeLLOs in that this 
source appears to have the lowest aperture mass of all VeLLOs. It 
also features the lowest N 2 H+ line width ([0.232± 0.048] km s" 1 ; 
the next higher value for VeLLOs is 0.3 kms -1 for L1521F). Fi- 
nally, LI 148 is the only "unevolved" VeLLO core in which inward 
motions have been detected. Since line asymmetries are easier de- 
tected for narrow lines (Eq. lB2t . this might be a consequence of the 
low N2H 4 " line width. 

5.4 Properties of VeLLO Outflows 

VeLLOs differ vastly in their outflow properties. For example, only 
2 out of the 7 VeLLOs discussed in Table|2](i.e., IRAM04191 and 
L673-7-IRS) drive a "prominent" extended outflow that can easily 
be detected using single-dish telescopes. The remaining VeLLOs 
all appear to drive outflows (based on Spitzer image morphology 
and interferometer observations), but sensitive high-resolution ob- 
servations are needed to detect these (like those presented in Sec- 
tion[3j. 

As seen in Table [2] the presence of prominent outflow activ- 
ity does not appear to correlate with any other dense core property. 
In the present sample, to be specific, an extended outflow is driven 
by a VeLLO in one of the most "evolved" cores, i.e. IRAM04191, 
which also has a large aperture mass. But also the VeLLO in L673- 
7 has a prominent outflow, despite a ve ry moderate aperture mass 
and failure in some of the ICrapsi et al] criteria. In addition, cores 
more evolved than IRAM04191 (e.g., L1521F), as well as cores 
very similar to L673-7 (e.g., L1014), do not drive such flow. This 
diversity could be a consequence of VeLLO evolution; VeLLOs 
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Table 2. VeLLOs and their Natal Cores 



Evolved Cores 



Unevolved Cores 



Property 


L1521F 


cna— MJVla 1 


1KAJV1U4 1 y t 


L673— 7 


L1014 


L328 


LI 148 


protostar properties: 
















outflow 


Y 


7 


Y 


Y 


Y 


? 


Y 


extended outflow 


11 


n 


Y 


Y 


11 


11 


n 




RJ 0.04 


[« 0.01] 


ss 0.07 


w 0.04 


ss 0.09 


» 0.05 


Kb 0.10 


Fi.scmd 2 /!®- 2 mjykpc 2 


0.3 


7 


0.3 


0.5 


1.9 


7 


sC 0.7 


dense core properties: 
















M4200 Au/A^0 


0.87 


[0.8] 


0.60 


0.39 


0.36 


[> 0.13/0.31] 


0.14 


(A r (H 2 )) 4 200 Au/10 22 x cm- 2 


3.0 


[2.7] 


2.0 


1.3 


1.2 


[> 0.4/1.1] 


0.5 


Af(N 2 D+) > 1.0 x 10 12 cm- 2 


Y 


[Y] 


7 


? 


11 


n 


n 


iV(N 2 H+) > 8.5 x 10 12 cm- 2 


Y 


[Y] 


Y 


■'! 


11 


Y 


n 


Af(N 2 D+)/A/(N 2 H+) > 0.1 


Y 


Y 


7 


? 


Y 


11 


7 


/ d (CO) > 10 


Y 


7 


? 


7 


n 


11 


n 


n c (H 2 ) > 2.5 X 10 5 cm" 3 


Y 


Y 


Y 


11 


11 


n 


n 


Au(N 2 H+) > 0.25 km s" 1 


Y 


Y 


Y 


Y 


Y 


Y 


n 


inward motions 


Y 


7 


Y 


n 


n 


ii 


Y 


r 70% < 4 800 AU 


Y 


Y 


Y 


Y 


Y 


Y 


Y 



;e et all Jl999Ll2009l);lAndre et alJJl999l);lMotte & Andrei l200ft;lBelloche et all <2002ll2006l):lYoung et alJj2004l) ; ICrapsi et alj 
5l l2006l) : lDunham et al l 12006. 20 1 Oah iKauffmann et alj j 2005ll2008L and this work): and lScaife et al] <2010l) 



Reference s. Data are fro m Lee i 
<2005allbl) : lBourke et all feOOl l 

Notes. The dense core properties are compared to the criteria for "evolved" cores bv lCrapsi etail l2005al) after conversion of the limiting density to our choice 
of dust emission properties. A 'Y' marker indicates that a criterion is met, a 'n' marker that it is not, a question mark that the data are insufficient to decide. 
Properties marked by square brackets ([...]) are not readily available from the relevant publications and must be gleaned from the other publis hed properties. 
For C ha-MMS 1 , the N 2 H + and N 2 D + column densities have to be derived from the abundances and dust column densities reported by Belloche et alj 
<2006t) . The aperture mass is derived from models for t he density structure presented in the same publication. For L328-IRS, the envelope masses are from 
lower limits of the 1.2 mm wavelength dust e mission jLee et alj|2009l ; 70 mjy within 2000 AU radius) and integration of the lLee et alj envelope density 
model to 4200 AU radius. We deviate from the lCrapsi et alj practice to infer inward motions from asymmetries in the N 2 H+ lines; instead, inward motions 
are either inferred from the analysis of line shifts (L1521F, L673-7, L1014, L328, and L1148), or from detailed modeling of line shapes (IRAM04191). For 
all sources the uncertainty in the internal luminosity is of order ±0.03 Lq. 



without outflows might be in a stage preceding or following an 
accretion phase, in which prominent outflows are produced. It is 
not clear whether other options — i.e., VeLLO formation without 
prominent outflows in any evolutionary stage — provide credible 
models. 

It is also not clear how the presence of prominent outflows re- 
lates to the concept of "episodic accretion", i.e., accretion (and thus 
outflow generation) occurring only duri ng short episodes of the 
protostellar lifetime dKenvon et alj|l990h . In this picture, VeLLOs 
would quasi-periodically create and gradually loose their promi- 
nent outflows. Dynamical ages for the outflows of IRAM04191 
and L673-7-IRS suggest that individual outflow periodjf] last 
a few 10 4 yr in VeLLOs. If VeLLOs without prominent out- 
flows are indeed in a dormant stage, this stage should prevail 
for an even longer timescale (> 10 4 yr), as indicated by the 
large fraction of VeL LOs without such outflows. Combination of 
Spitzer source counts fevans et al.l2009l) with evolutionary models 
jDunham et al.l2010bl) indeed shows that YSOs might accrete most 
of their mass in < 10 bursts occurring over a period ~ 2 X 10 yr. 
These models do thus imply dormant periods of several 10 4 yr, just 
as suggested by VeLLO outflow statistics. 



6 DISCUSSION: L1148-IRS IN CONTEXT 
6.1 Evolutionary Stage & Mass of L1148-IRS 

L1148-IRS has a bolometric temperature (110 ± 15) K (Section 
I4.lt . This implies that LI 1 148-IRS bel ongs to the protoste llar class 
I, and suggests an age < 7 x 10 5 yr jEvans et alj|2009l) . We use 
this limit in the following. 

As shown by lYoung et alj | |2004|) , the internal luminosity of ^ 
0.1 Lq limits the present mass of all VeLLOs to a most likely value 
< 0.1 A/q, if the luminosity is dominated by accretion power. This 
is probably the case for LI 148-IRS, since it appears to drive an 
(supposedly accretion-powered) outflow. Further, the young brown 
dwarf F U Tau (a class II obj ect of < 10 6 yr) has a luminosity of 
0.2 L dLuhman et al]|2009l : > 0.1 L & for the brightest compo- 
nent, if FU Tau is an unresolved binary ). This object might be un- 
usually luminous dLuhman et alj 120091) . but its properties support 
the notion that YSO luminosities < 0.1 Lq imply a mass in the 
brown dwarf domain. Thus, we find that the most likely present- 
day mass is < 0.1 Mq. 

Comparison with evo lutionary tracks bv lBaraffe et alj (l998) 
and IChabrier et alj ( feood) suggest a very conservative mass limit 
■ C 0.35 (when r educing the model luminosities by a factor 
3: lReiners etal]|2005l) . Models bv lWuchterl & Tscharnuten d2003h 
suggest < 0.1 Mq, if one only considers early stages with sig- 
nificant accretion. Unfortunately, their models are not well tested 
against observations. 



5 Prominent outflows are detectable during this period. This does, however, 
not require that they are continuously driven during this time. 



To understand the collapse conditions in LI 148, consider the col- 
lapse of a hydrostatic isothermal sphere of radius R that is subject 



12 J. Kauffmann et al. 



to an environmental pressure, P en v To be unstable, the sphere must 
have a mass 



M > 2.0 M 







R 



WKJ VO-ipc 
and the pressure must obey 

Pxt 1 n 4 —3 I -^s 

— — > 3.3 x 10 cm 



r 



(i) 



(2) 



,iok/ VO-ipc 

jBormojl 19561 and lEbert|[l95i k B is Boltzmann's constant). Star 
formation requires that the inner parts of the protostellar en- 
velope are unstabl^f]. We can use the aforementioned aperture 
mass, M4200AU, to test this condition. For R = 4200 AU and 
T g = 10 K, Eq. G} yields the condition M > 0.41 Mq. Thus, 
M4200AU fails the condition for gravitational instability. The reser- 
voir actively forming the VeLLO must therefore have a radius 

< 4200 AU. This implies an unstable mass smaller M4200AU, i.e. 

< 0.14 Mq (Section|2T}. Radii < 4200 AU also imply a pressure 
Pext/ks > 8 x 10 5 cm -2 . As we discuss in Section [63l this pres- 
sure is relatively high. This has implications on VeLLO formation 
models. 

Combined with the present-day mass of L1148-IRS, < 
0.1 Mq, this estimate of the collapsing mass suggests a final mass 

< 0.24 Mq. If this VeLLO has a present-day mass < 0.05 Mq, 
and the immediate envelope of < 0.14 M q is accreted with an ef- 
ficiency ~ 1/3 (as deemed typical; e.g.. lAlves et alj|2007l) . then 
L1148-IRS could have a final mass ev en below the brown dwarf 
limit. After L328-IRS dLee et alj|2009l) . this is the second VeLLO 
for which such a small final mass is derived. 



6.2 The Nature and Evolution of VeLLOs 

There is probably general consent that VeLLOs have a most likely 
present-day mass < 0.1 Mq (Section [64}. The final mass at the 
end of the accretion phase might, however, be much larger. 

VeLLO natal cores differ significantly in the mass of gas at 
high density they contain (Section[5]l. Since VeLLOs form from this 
gas, these differences ma y suggest that differ ent VeLLOs take dif- 
ferent evolutionary paths dBourke et al.l2006h . Some may have sig- 
nifica nt future accretion and attain a stellar final mass l l Andre et al.l 
1 19991) . while lack of accretion in others could render these ob- 
jects young brown dwarfs dYoung et"allll2004l ; iHuard et all 1 120061) . 
Table [2] thus presents a first attempt to systematically characterize 
the likeliness of these different evolutionary paths for several VeL- 
LOs. The number of positive indicators for significant amounts of 
high density gas near the VeLLO increase towards the left of Ta- 
ble^ The VeLLO final mass is therefore likely to increase in this 
direction too. L328-IRS and L1148-IRS hold a peculiar position 
among VeLLOs, given that their most likely final mass is very low 
« 0.15 Mq and < 0.24 Mq, respectively; see iLee et al]|2009l 
and SectionHTl- 



6.3 Formation of Brown Dwarfs: Lessons from L1148 

One unsolved aspect of isolated brown dwarf (BD) formation is the 
large difference between the final BD mass, Mbd < 0.08 Mq, and 

B To be specific, this is a requirement only for models with inside-out col- 
lapse. The evolut ion of these collapsing spheres is simi lar to those stud- 
ied by [Shu| {l977l). One major difference, though, is that lShul considers — 
strictly speaking — infinite spheres with singular density profiles. His struc- 
tures are unstable to perturbations at all radii, just as required here. 



the expected mass of the collapsing natal dense core, M: in most 
cases M Mbd- L1148-IRS can be used to study this problem, 
even in case this VeLLO has a final mass exceeding the BD mass 
limit. 

Consider once again pressure-bound hydrostatic isothermal 
spheres. Combination of Eqs. (Q~|[2]l yields the condition 



M > 1.2 Mr. 



© 



10 K 



-Penv / 

10 5 Kcm- 



-1/2 



(3) 



M > 1 Mq . Ty pical values of r\ are believed to be ~ 1/3 (e.g.. 

stellar masses ~ 0.3 Mq. 



for unstable spheres. This system will form a star (or BD) of mass 
M* = r)M, where 77 ^ 1 is the star formation efficiency. For 
typi cal temperatures (~ 10 K ) and pressures (a few 10 5 Kcm -3 ; 
e.g. iBertoldi & McKeeill992l) , one obtains the collapse condition 
M > IMq. Ty pical values of rj i 
lAlves et al.ll2007l) , yielding typical : 

In this framework, to achieve lower stellar masses, BD for- 
mation requires unusually small temperature (T g -C 10 K), large 
pressure (Penv /fas > 10 s Kcm -3 ), low efficiency (77 < 1/3), 
or a combination of these. Ll 148-IRS provides a test case for such 
reasoning: Eq. Q} suggests M < 0.14 Mq for the collapsing core 
in this region. 

The structure of the Ll 148 core thus appears to require pecu- 
liar environmental conditions. We have no satisfying explanation 
why and in which sense (temperature, pressure, etc.) these should 
prevail in L1148. One can, however, speculate that the unusual 
kinematic conditions in this region might have played a role: mul- 
tiple velocity components (in C ls O), velocity jumps near Ll 148- 
IRS (in CCS) and contraction motions on spatial scales > 0.3 pc 
are not commonly observed in dense cores (see Section [2~T1 and 
Appendix |B3j- 

The en vironmental condit ions in L1148 have actually been 
explored bv lNutteretal] d2009h . Spatial offsets between intensity 
peaks at different wavelengths are consistent with external heat- 
ing of the Ll 148 core. On e could imagine th at the core is sculpted 
by such a radiation field. iKirk et all (2009) suggests indeed that 
the entire Cepheus Flare Complex is shaped by external triggering. 
Related externally triggered compression events would be suited to 
produce low critical masses via an increase in P on v 

It might, of course, also be that the initial unstable mass of the 
dense core was ~ 1 Mq, as suggested by Eq. (|3j, but that we now 
observe a much more evolved p hase, where the initial mass is no 
longer relevant. ILee et all I d2009h . for instance, suggest that L328- 
IRS might have moved out of its natal core. This wou ld leave this 
VeLLO with a small mass reservoir to accrete from. IHuard et"al] 
d2006h present a similar argument for L1014-IRS. The complex 
kinematics in Ll 148 might suggest such a scenario also for Ll 148- 
IRS. In fact, Section |2~T| shows that there is a possible — though not 
significant — 7'.'6 spatial offset (i.e., 2500 AU) between Ll 148- 
IRS and the core's column density peak. The growth of Ll 148-IRS 
might thus be stinted, because it has detached from the dense core 
from which it cannot accrete from it any longer. 

Similarly, Ll 148-IRS could be in a very advanced stage of 
evolution where most of the envelope has been dispersed. In prin- 
ciple, strong VeL LO outflows (e.g., as those from IRAM04191; 
lAndre et al ll 19991) might be able to do this. 

The credibility of these scenarios is unclear, though. Consider 
L1148-IRS as an example. Section [6T| suggests that the present- 
day mass of this VeLLO is < 0.1 Mq. The immediate environment 
of Ll 148-IRS contains ~ 0.14 M (from the 4200 AU aperture 
mass). If outflows have dispersed the gas, they must thus have dis- 
rupted a core of about 0.75 Mq. It is not known whether dispersal 
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by outflows is efficient enough to do this dArce et alJ,l2007t) . Sim- 
ilarly, it is not obvious how VeLLOs should be able to leave their 
parental dense core. That would require that accretion from the gas 
reservoir onto the VeLLO increases the relative velocity between 
the gas reservoir and the VeLLO. This seems unlikely. 



7 SUMMARY 

We present a detailed analysis of a candidate VeLLO in the 
Cepheus Flare, LI 148-IR S. This work goes beyon d the initial an- 
nouncement of discovery (Kauffmann et al. 2005;|) by presenting: 
(i) a comprehensive and sensitive search for a CO outflow; (it) de- 
tailed modelling of the observed spectral energy distribution; and 
(Hi) extensive discussions of the properties of the natal dense core. 
In addition, (iv) the analysis is augmented by presenting LI 148- 
IRS in the context of an up-to-date VeLLO catalogue. The main 
conclusions are as follows. 

(1) LI 148-IRS drives a compact (~ 5", about 1800 AU in pro- 
jection) CO outflow (Section [3}. This confirms LI 148-IRS as a 
clear-cut protostar, and associates it with LI 148 through its veloc- 
ity. The outflow can only be detected in interferometer maps, such 
as those from the Plateau de Bure Interferometer presented here; 
single-dish data, e.g. from the IRAM 30m-telescope, do not reveal 
it. Weak compact (~ 2") probable outflow nebulosity is seen near 
2(im wavelength. 

(2) Modelling of the spectral energy distribution suggests the 
presence of a point source of 0.08 to 0.13 Lq and T c g = 
2500 to 5000 K effective photospheric temperature (Section |4](. 
Strong emission near 100u.m wavelength requires this source to be 
deeply embedded in the natal dense core; it corresponds to a source 
of infrared class or I. LI 148-IRS does therefore meet all criteria 
to be a VeLLO (i.e., L ^ 0.1 Lq and embeddedness in a core), 
though this is just marginally the case for the luminosity. 

(3) The natal dense core is characterized by densities and column 
densities that are unusually low for a star-forming core (Section 
12 . lb - Multiple velocity components, velocity jumps, and velocity 
gradients are detected in the core, pointing towards unusual kine- 
matic conditions in this region. Strong evidence for contraction mo- 
tions is found in asymmetries of optically thick lines. 

(4) Among all VeLLOs known, LI 148-IRS stands out as the one 
with the lowest envelope mass and a core of unusually low (col- 
umn) density (Table[2]m Section[6](. It is a riddle how a star could 
have formed under these conditions. 

As in all other VeLLOs, the low luminosity suggests a substellar 
present mass for LI 148-IRS, i.e. < 0.1 Mq. The gravitationally 
unstable mass surrounding LI 148-IRS appears to be < 0.14 Mq. 
Their sum suggests a very low final mass, so that LI 148-IRS may 
be a good candidate for a proto brown dwarf forming in isolation. 
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APPENDIX A: PHOTOMETRY FOR L1148-IRS 

Table |Al] presents the SED of LI 148-IRS, compiled from data 
sources listed in Section |4~T1 Here we give some further process- 
ing details. 

In the IRAC and the MIPS 1 band, we average results, if more 
than one flux density measurement is available. For the MIPS2 
and MIPS3 band, we use our own aperture photometry results in- 
stead of the values derived by point-source fitting; this is likely to 
be better suited for the slightly extended emission. The SHARC- 
II and SCUBA upper limits are derived by adding the observed 
flux and the no ise level times t hree (our SHARC-II limits super- 
sede those from lWu et alj|2007l). We chec k our photometry against 
further Spitzer M IPS jKirk et all l2009h and independent Akari 
jNutter et alj2009t) studies of the regior|J. The data are usually con- 
sistent within their uncertainties. A notable exception is the Akari 



flux of 130 mjy at 90 fim wavelength; this is inconsistent with the 
much higher MIPS flux densities at 70 and 160u.m. We have no 
explanation for this, but elect to rely on our higher quality Spitzer 
data. 



APPENDIX B: THE L1148 DENSE CORE 

As already explained in Section 12.11 iKauffmann eTai] l l2005h 
lacked the space to describe their data in detail. To provide a proper 
reference, we do therefore present those observations in this ap- 
pendix. We also discuss new observations of contraction motions 
in LI 148 (Appendix |B5t . To improve the readability of this paper, 
their rather straightforward discussion is removed from the main 
text. 

In this appendix, we do therefore first introduce a previously 
rather undocumented data set that had never been documented well, 

1. e., our Effelsberg observations (Appendix IB lb . We then discuss 
the spatial mass distribution of the core in more detail than appro- 
priate in the main text (Appendix lB2l >. This is foll owed by extended 
discus sions of the line emission data first used by I Kauffmann et all 
fcOOSl) (Appendix [B3] velocity field of the cloud; Appendix IB4I 
molecular abundances). We finish with an analysis of new data on 
contraction motions in the cloud (Appendix |B5| >, 

Bl Additional Data: Effelsberg lOOm-telescope 

Spectroscopic observations with the Effelsberg lOOm-telescope 
(operated by the Max-Planck-Institut fur Radioastronomie in 
Bonn, Germany) were done 2005 April 13 and April 29 to May 

2. The (J, F — 2, 1-1, 0) transition of CCS was observed with the 
facility 1.3 cm primary focus spectral line receiver and the AK90 
autocorrelator as the backend. All data was taken in the frequency 
switching mode. The pointing corrections in azimuth and eleva- 
tion were mostly below 6", but corrections of 15" occurred once. 
The focus corrections were below 0.6 mm. Pointing was checked 
at least once in 4 hours and the focus corrected at least once in 6 
hours. Further details are listed in Table IB 1 1 

The spectroscopic observations were reduced using CLASS. 
We removed low-order baselines and partially used Fourier (FFT) 
filtering to further improve the baseline quality. Further details of 
the data reduction are given in the respective sections below. 

The data were calibrated using cross-scans towards 3C48, 
3C123, 3C147, 3C286, and NGC7027. The observed telescope 
respon se is combined with the flux densities given by lOtt et al.l 
d 19941) to determine the calibration factor. We use a factor of 
3.62 to convert the dimensionless telescope counts into main 
beam brightness temperature. Please refer to \protect\vrule 
width Opt \protect \href { http : / / www . mpif r-bonn . 
mpg . cle/ staff /tpillai/ef f_calib/ ef f_calib . 
html } { http : / / www . mpif r-bonn . mpg .cle/ staff/ 
tpillai/ef f_calib/eff_oalib . html } for more details. 

B2 Spatial Mass Distribution 

Based on the dust em ission maps shown in Figures [TJ and [2] 
IKauffmann et alj J2008I) noted that the LI 148 region contains two 
parallel filaments of 3 Mq and 9 Mq, respectively of 1' x 5' size 



7 D. Nutter provided us with dedicated photometry, based on the data from 



his Akari study jNutter et al. 2009|). We are grateful to him and J. Kirk for 
providing us their results in advance of publication. 



16 J. Kauffmann et al. 

Table Al. The spectral energy distribution of LI 148-IRS 



Band Wave Length Bandwidth Resolution Apertures Flux Density Calibration Uncertainty 

(im (im arcsec arcsec mjy 
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Notes. For every band the central wavelength and width of the filter, the telescope resolution, the apertures used for photometry, the measured flux density 
including color corrections (where necessary) and its uncertainty due to noise, and the calibration uncertainty are given. For most bands photometry is done by 
point-source fitting. Then the aperture to that the flux density refers is given by the point spread function (PSF) and we mark these bands by listing "PSF" as 
the used aperture. For the other bands, where we perform aperture photometry, the first number given for the aperture is the diameter across which the intensity 
is integrated, while the given range indicates the diameter of the inner and outer circle enclosing the area used to obtain an estimate for the background intensity 
(does not apply to bolometer data). The uncertainty in the MIPS 3 band reflects the flux density variation due to different aperture choices. 



Table Bl. Spectroscopic observations at the IRAM 30m-telescope and the Effelsberg lOOm-telescope 
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J = 3-2 
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12 CO 


J = 2-1 
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10 
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377 f 
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c 18 o 


J = 2-1 
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HCO+ 


J = 1-0 
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0.64 ± 0.04 
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H 13 CO+ 


J = 1-0 "primary" 
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0.20 ±0.04 
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HCN 


J,F=1, 1-0, 1 
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J, F = 1, 2-0, 1 
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0.16 ±0.02 
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J, F = 1, 0-0, 1 


88.6339360 1 


0.09 ±0.02 






. . ditto 







References. (1) iMuller et al l 120051) ; (2) bore et alj 120041) : (3) iGottlieb et all J2003ri ; (4) Iwinnewisser et alj dl997h ; (5) ICazzoli et al] J2003l) : (6) 
ISchmid-Burgk et al l <2004h 

Notes. a peak main beam brightness temperature per spectrum (and noise level measured as the root mean square deviation at velocities free of line emission) at 
the position closest to LI 148-IRS, determined using Gaussian fits to the lines, where dashes ("-") indicate non-detections, and upper limits are the noise level 
times three; b uncertainty in frequency, and implied uncertainty in velocity; c frequency throw in frequency-switched observations; d system temperature, in 
main beam brightness scale; c velocity resolution of the spectral observations, after smoothing if applicable; f for mapping observations in December 2004; s 
for pointed observations in July 2005 
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(i.e, 0.1 x 0.5 pc 2 ; both filaments have a similar size). The north- 
western filament contains LI 148-IRS and has the lower mass. The 
H2 peak column densities are 6 x 10 cm and 8 x 10 cm 

To derive radial density laws, we execute a new analysis of 
cuts through the dust emission peak near LI 148-IRS. They are 
taken at a position angle 135 deg, i.e., perpendicular to the fila- 
ment axis. The cuts (i.e., dust emission intensity as a function of 
offset from the peak) are compared to intensity profiles expected 
for density laws g(r) oc -a and constant dust temperature. The syn- 
thetic intensities are calculated taking the limited telescope resolu- 
tion into account. Within the uncertainties, the dust emission distri- 
bution is consistent with pure power-laws with a — 1.75 to 2.0. 
In other words, compared to the beam diameter, regions of constant 
density (i.e., a = 0) cannot have a significant size. 



B3 Velocity Field 

iKauffmann et al.1 J2005I) used CCS (J F = 2i-l ) maps (not re- 
peated here) to gauge the velocity field of the LI 148 region on very 
large spatial scales (i.e., 3> 1', or 3> 0.1 pc). The maps reveal 
similar velocities ~ 2.7 kms - for the two dust filaments (that 
are traced by the CCS intensity), and full width at half maximum 
(FWHM) line widths (from Gaussian fits) of 0.15 to 0.25 kms" 1 . 
There is a velocity gradient along the north-western dust fil- 
ament harboring LI 148-IRS: the filament's north-eastern edge 
(that contains LI 148-IRS) has a velocity 2.63 kms" 1 , while the 
south-western end moves with 2.78 kms -1 . This difference of 

0. 15 km s - is similar to the line width, and thus rather significant. 
The transition from "low" to "high" velocity occurs rather abruptly, 

1. e., over a distance of only « 1', or ~ 0.1 pc. This suggests that 
this filament actually consists of two physically independent ob- 
jects that form a single object only in projection. The line width is 
slightly elevated in the transition region, but not to a level which 
would clearly indicate the physical interaction of the two velocity 
components (e.g., not by a multiple of the average line width). The 
transition zone lies only ~ 1' to the south-east of LI 148-IRS, and 
is thus relatively cl ose to our target object. 

In the maps of IKauffmann et al] <2005h . smaller spatial scales 
are well probed by the C ls O (J = 2-1) maps of size 100" x 100", 
or (0.16 pc) 2 . As shown in Figure IbTI they reveal two velocity 
components in the immediate vicinity of LI 148-IRS: one at veloc- 
ities ~ 2.6 kms" 1 , and one at ~ 3.0 kms" 1 . The FWHM line 
widths (from Gaussian fits) show significant variation across the 
map; they are 0.2 to 0.5 kms" 1 , when only taking reliable fits into 
account. The velocity of the lower-velocity emission globally in- 
creases from north-east to the south-west, but this trend is not as 
distinct as seen in CCS. Interestingly, the C 18 integrated intensity 
presented in Figure lB2l does not correlate with the dust intensity 
distribution (also see Appendix IB4t . 

The kinematic s of the very dense gas a re traced by the N2H + 
(J = 1-0) maps oflKauffmann et al. (2005). As seen in Figure lB2l 
these cover the same region as the C 18 maps. The spectra show no 
signs of multiple velocity components. Characteristic velocities can 
be derived by simultaneously fitting all N2H" 1 " hyperfine satellites. 
They are of order 2.6 kms -1 . Velocities increase from north-east 
(~ 2.55 kms" 1 ) to south-west (~ 2.63 kms" 1 ), just as seen for 
CCS and C 18 0. The FWHM line widths are 0.15 to 0.25 kms -1 , 
with a peak towards LI 148-IRS. 

The low observed line widths < 0.25 kms -1 are notewor- 
thy; N 2H + lines in supposedly "evolved" cores l lCrapsi et al.l 
l2005al) . as well as CCS transitions in star-forming cores 
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Figure Bl. Spectra in the C ls O J = 2-1 transition towards LI 148-IRS as 
observed with frequency switching on a 5 x 5 point grid with 30" spacing. 
Offsets are given relative to the approximate dust emission peak position 
(a = 20 h 40 m 56?5, S = +67° 22' 58'.' [J2000.0]). The LSR velocity 
and intensity in the T m ^ -scale are shown for the spectrum in the lower left 
corner. Gaussian fits to velocity components are drawn for those positions 
where the component identification is obvious. Small numbers give for ev- 
ery component fitted the derived LSR velocities (upper lines) and line width 
(lower lines; both in km s — 1 ). 



llde Gregorio-Monsalvo et alEoOfih . are observed to typically have 
lines wider than 0.25 kms" 1 . It is, e.g., possible that the wide lines 
in non- VeLLO cores are a manifestation of large-scale rapid in- 
ward motions associated with accretion on forming stars. If these 
motions are slower in VeLLO cores, or their spatial extent is 
smaller, then if would be understandable why VeLLO cores have 
so narrow lines. However, in the light of the present data situation, 
any such consideration remains speculative. 



B4 Molecular Abundances 

The CCS integrated i ntensity does not correla te well with the dust 
emission (Figure 3 of lKauffmann et alj|2005l) . In contrast, N 2 H+ 
(1-0) emission (Figure lB2l traces the core very well. 

The correlation between the C ls O (2-1) emission and the dust 
distribution is unusually weak (Figure lB2t : while the dust emission 
filament runs NE-SW,the C 18 emission features a distinct wedge 
that starts near LI 148-IRS and flares towards SSE (position angle 
roughly 160°). This structure is part of the aforementioned lower- 
velocity C 18 component: it is still seen if subtracting the emission 
from the higher-velocity component from the map (by using the fits 
shown in Figure iBll . 

Kauffmann et al. I d2005h calculate abundan ces using the proce - 
dures and molecular constants laid out by ICrapsi et alj d2005al) : 
for N2H + , optical depths and excitation temperatures are de- 
rived from hyperfine-fits, while C ls O lines are assumed to be 
optically thin and to have an excitation temperature of 10 K. 
Under these assumptions, the CO depletion factor (which is 
supposed to increase when CO freezes out onto dust), /<j = 
[A r (C ls O)/A r (H 2 )]o/[A r (C 18 0)/A''(H2)]obB, is calculated from 
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Figure B2. Maps of C ls O (J = 2-1) and N 2 H+ (J = 1-0) integrated 
intensity. Grid markers give the positions observed in our study. Shading 
illustrates the intensities recorded during the observations. Contours out- 
line the dust emission detected by MAMBO. For N2H+, marker sizes and 
numbers indicate the radial velocity derived from hyperfine fitting; crosses 
are drawn if the emission is too faint to be fitted reliably. A star marks the 
position of L1148-IRS. 
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Figure B3. Spectra towards the positions (left to right) and transitions (top 
to bottom; J, F = 1, 2-0, 1 for HCN, J, F\,F = 1, 0, 1-0, 1, 2 for 
N2H + ) defined in Table lB2l Dotted vertical lines indicate the systemic 
velocity derived from the N 2 H+ (J, Fi,F = 1,0, 1-0, 1, 2) transition. 
All well detected lines believed to be good tracers of infall have peaks 
blueshifted from the systemic velocity, indicating inward motions. 

The CO depletion factor peaks towards L1148-IRS, where 
/d = 7.8 ± 1.3. Compared to other dense cores, this is relatively 
low, in particular since this approach is likely to overestimate the 
true value (see above). For the same position, an N2H + column 
density of (3.1 ± 0.6) x 10 12 cm -2 is derived from deep new in- 
tegrations to probe contraction motions (Appendix |B5t . We do not 
detect significant N 2 D + emission. Assuming that the lines are op- 
tically thin, we derive N 2 D + column densities < 0.5 x 10 cm~ 2 
(using the intensity noise level times three). 



the observed (subscript 'obs') and "canonic" ('0') dust and CO col- 
umn densities, N, as 



fa = 0.42 



Sbeam 
1.2 mm 



(mj 



jr mb (c 18 o[2-i])d« 

Kkms- 1 



,(B1) 



where Si.^Sm is the dust emission intensity (per 11" beam), and 
/ r m b(C ls O[2 - 1]) dy is the velocity- i ntegra ted C 18 intensity 
in Tmb-scale. Note that lKauffmann et aU l l2005t) use — in deviation 
to all other parts of this paper — dust emiss ion intensity-to-mass 
conversio n factors fromlCrapsi et alj instead of those sug- 

gested bv lKauffmann et alj ( 120081) . This scheme is adopted to ma ke 
our depletion factor s more comparab le to those of ICrapsi et all If 
one were to use the iKauffmann et au recommendations, /a would 
decrease by a factor 2.54. Comparison shows that this approach 
typically overestimates the lCrapsi etafl J2005d ; their Table 6) esti- 
mates from (1-0) line data by a factor of 1.3, probably because of 
different excitation conditions for the (1-0) and (2-1) lines and dif- 
ferences in beam matching (for the present data, the beam in C ls O 
observations is smaller than the one in the smoothed MAMBO 
maps used to get Si.fmm)- 



B5 Contraction Motions 

The systemic reference velocity of the dense core, v rc { , can be well 
measured using the optically thin J,F\,F — 1, 0, 1-0, 1, 2 tran- 
sition of the N2H 4 " dense gas tracer. Towards L1148-IRS, we de- 
rive a velocity of 2.60 kins" 1 (Table lB2t . Many supposedly opti- 
cally thick lines do, however, peak at lower velocities (Figure lB3l . 
i.e. they are blue-shifted. Such veloc ity shifts are characteri stic of 
dense cores undergoing contraction dLeung & Brownlll977t) . The 
velocity offset of a line peaking at velocity Vi can be characterized 
using dimensionless velocity differences, 



Svi 



V le { 



Aw re f 



(B2) 



dMardones et alj[l997l) . where At; is the line width at half peak 
intensity. If Svi < 0, then the line is blue-shifted, which suggests 
contraction motions. Expansion implies Svi > 0. 

Table IB2l gives observed peak velocities in LI 148. These are 
derived by fitting Gaussian lines to the peaks (and sometimes ad- 
ditional Gaussians to the base). Table |B3l gives the corresponding 
dimensionless velocity differences. As outlined in Figure [2] we 
searched three positions for contraction motions, in order to remove 
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Table B2. Line Velocities and Widths for several Transitions observed towards selected Positions. 



Line 


LI 148 


-IRS 


offset position 


southern peak 




20 h 40 m 56?5, H 


-67 o 22'58'/0 


20 h 41 m 01?7, ±67°23'28'.'0 


20 h 41 m 14?0, ±67°20'31('5 




fLSR/km s" 1 


Aii/km s _1 


DLSR/kms 1 


Aii/km s 1 


"LSR/km s 1 


Ati/kms 1 


N 2 H+ (J, Fi,F = 1,0, 1-0, 1, 2) 


2.602 ±0.015 


0.244 ± 0.030 


2.583 ± 0.023 


0.232 ±0.048 


2.614 ±0.013 


0.243 ± 0.029 


CS (J = 2-1) 


2.418 ±0.005 


0.149 ± 0.017 


2.488 ±0.003 


0.118 ±0.010 


2.509 ±0.002 


0.119 ± 0.004 


CS (J = 3-2) 










2.513 ±0.005 


0.106 ± 0.015 


HCO+ (J = 1-0) 


2.218 ±0.001 


0.204 ±0.013 


2.327 ±0.003 


0.119 ±0.008 


2.356 ±0.002 


0.115 ± 0.006 


H 13 CO+ (J = 1-0) 


2.681 ±0.016 


0.420 ± 0.038 


2.652 ±0.015 


0.314 ±0.028 


2.628 ± 0.008 


0.293 ± 0.017 


HCN(J,F = 1,2-0, 1) 


2.360 ±0.013 


0.197 ±0.048 


2.398 ± 0.006 


0.127 ±0.018 


2.424 ± 0.008 


0.125 ± 0.021 


HCN (J, F = 1, 1-0, 1) 






2.429 ±0.018 


0.198 ±0.046 


2.438 ±0.012 


0.113 ± 0.032 


HCN(J,F = 1,0-0, 1) 


2.601 ±0.041 


0.566 ± 0.086 


2.460 ±0.016 


0.153 ±0.033 


2.469 ±0.010 


0.124 ±0.034 



Notes. The observations hold for the J2000.0 positions defined in the header. These are also marked in Figure [2] If more than one line has to be fitted 
simultaneously to achieve a good fit to the narrow peak and the broad base of a line, then the parameters for the component with the lower velocity are listed. 
No numbers are given for transitions not detected well enough to allow fitting. The lines believed to be optically thick all peak at velocities below the systemic 
one derived from N2H+. 



Table B3. Dimensionless Velocity Shifts towards selected Positions 



Line 



CS 

J = 2-1 
J = 3-2 

HCO+ 
J = 1-0 

H 13 CO+ 
J = 1-0 

HCN 

J, F = 1, 2-0, 1 
J, F = 1, 1-0, 1 
J,F = 1,0-0,1 



SVi = (Vi - f N2 H+)/ At, N 2 H+ 
LI 148-IRS offset position southern peak 

-0.75 ±0.09 -0.41 ±0.08 -0.43 ± 0.05 

-0.42 ±0.05 



-1.57 ±0.19 -1.10 ± 0.23 -1.06 ±0.13 



0.32 ±0.04 0.30 ±0.06 0.06 ± 0.01 



-0.99 ± 0.12 



0.00 ±0.10 



-0.80 ±0.16 
-0.66 ± 0.14 
-0.53 ± 0.11 



-0.78 ±0.09 
-0.72 ±0.09 
-0.60 ±0.07 



Notes. Dimensionless velocity shifts in several transitions towards the po- 
sitions defined in Table IB2I No numbers are given for transitions not well 
enough detected to allow fitting. The lines believed to be optically thick all 
do have negative line shifts indicative of inward motions. 



depth (in LTE, the optical depth ratios for different F are 1:5:3 for 
[0—1] : [2— 1 ] : [ 1— 1]), just as expected for contracting cores. Towards 
LI 148-IRS, the possibly optically thin F = 0-1 line of HCN has 
an unusually large line width of 0.57 kms -1 . This may hint at me- 
chanical energy input by the YSO, for example in the form of out- 
flows. The H 13 CO + (1-0) is supposedly optically thin and peaks 
at higher velocities then the optically thick lines. This confirms the 
velocity differences with respect to N2H + , although the line cannot 
serve as a precise velocity reference becaus e of unresolved hyper- 
fine substructure dSchmid-Burgk et alj|2004) . 

It thus appears that the natal dense core of LI 148-IRS is in- 
deed contracting. This is suggested by the negative velocity differ- 
ences and the particularly strong line asymmetries towards the em- 
bedded source, LI 148-IRS. The prevalence of blue-shifts through- 
out the LI 148 region may actually indicate that the entire region 
(i.e., on scales larger than the ~ 0.3 pc containing the three posi- 
tions probed) is contracting. 



biases from possible foregrounds to LI 148. All lines are suppos- 
edly optically thick. This only excludes the "isolated" N2H + hy- 
perfine satellite (i.e., J,Fi,F = 1, 0, 1-0, 1, 2) and the H 13 CO+ 
(1-0) transition, which can therefore work as velocity references. 

All well-detected optically thick lines peak at velocities blue- 
shifted from the N2H + reference velocity. This only excludes the 
(F = 0-1) emission of HCN towards LI 148-IRS, which is a 
particular case discussed below. In all tracers, the dimensionless 
line asymmetries are largest towards LI 148-IRS. This position is a 
"strong infall candi date", following the criteria for CS (2-1) emis- 
sion put forward bv lLee et alj J 19991 : i.e., a CS emission peak with 
significant blue-shift; red-shifted emission at a brightness below 
the one of the blue-shifted emission; and no asymmetries in the 
isolated N 2 H + line). Also, the values of the CS and HCN line 
asymmetries towards LI 148-IRS are unusu ally large for sta rless 
cores: only 18% of the CS cores probed by iLee et al.1 dl999h . re- 
spectivel y 40% of the HCN cores studied in the (F = 2-1) tran- 
sition bv lSohn et al. I d2004) . have l arger asymmetri es (for CS, we 
use the data after the reanalysis bv lLee et alj|200ll) . Note that the 
asymmetry in the HCN satellites decreases with decreasing optical 



